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Abstract 

Astronomical observations have shown that small carbonaceous molecules can 
persist in interstellar clouds exposed to intense ultraviolet radiation. Current astro- 
chemical models lack quantitative information on photodissociation rates in order 
to interpret these data. We here present ab initio multi-reference configuration- 
interaction calculations of the vertical excitation energies, transition dipole moments 
and oscillator strengths for a number of astrophysically relevant molecules: C3, C4, 
C 2 H, I- and c-C 3 H, I- and c-C 3 H 2 , HC 3 H, Z-C4H and Z-C 5 H. Highly excited 
states up to the 9'th root of each symmetry are computed, and several new states 
with large oscillator strengths are found below the ionization potentials. These data 
are used to calculate upper limits on photodissociation rates in the unattenuated 
interstellar radiation field by assuming that all absorptions above the dissociation 
limit lead to dissociation. 
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1 Introduction 



Of the more than 130 different molecules found in interstellar space, an impor- 
tant class is formed by the unsaturated carbonaceous species. Carbohydrides 
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ranging from small molecules such as C2H (Q]) to long chains like CsH (0) and 
HCnN (jil) have been detected through their millimeter transitions in cold dark 
clouds like TMC-1 for several decades. Some of these molecules have also been 
seen in their cyclic form, with C-C3H2 as the best-known example (0). Even 
negative ions, in particular CeH~, have now been detected (ja). Bare carbon 
chains are likely present as well but do not have a permanent dipole moment 
and can therefore not be observed through their pure rotational transitions in 
emission. Instead, the smallest members of this family, C2 and C3, have been 
detected in diffuse interstellar clouds, -i.e., clouds which are not completely 
opaque to visible and ultraviolet radiation- through their electronic absorp- 
tions against bright background stars (Ji!; 0; H) . c— C3H2 has been seen in diffuse 
clouds as well, through absorption at radio wavelengths against background 
quasars (Q). The relatively large abundances of these non-saturated molecules, 
in spite of the fact that there is 10 4 times more hydrogen than carbon in inter- 
stellar clouds, are a vivid demonstration that interstellar chemistry is not in 
thermodynamic equilibrium. Instead, the kinetics of the reactions that form 
and destroy these molecules need to be taken into account explicitly in order 
to explain their abundances. In cold dark clouds, these traditionally involve a 



series of ion-molecule and neutral-neutral reactions (lldl ; llll ) 



More recently, several of the carbohydrides C n H m have also been observed in 
so-called Photon-Dominated Regions (PDRs) (Il2l ). i.e., clouds which are ex- 
posed to intense ultraviolet (UV) radiation. Species like C2H and C-C3H2 have 
been detected in the Orion Bar PDR (1131 ) and even more complex molecules 



like C3H2 and C4H have been seen in the Horsehead nebula, M38 (114j ; Il5l ). 
Together with the above mentioned observations of diffuse clouds, the data 
demonstrate that these molecules can exist in UV-exposed environments with 
abundances comparable to those in cold dark clouds. Traditional PDR chemi- 
cal models cannot explain the high abundances of these molecules, leading to 
speculations that they are perhaps produced by fragmentation of even larger 
carbonaceous molecules such as polycyclic aromatic hydrocarbons (PAHs) 
(flol). However, a major uncertainty in these and other models are the pho- 
todissociation rates of the carbohydrides, which are basically unknown. 

Other regions in which ultraviolet photons play a role in the chemistry in- 
clude the envelopes of evolved stars and cometary atmospheres. The same 
carbohydride molecules have been detected in the envelope of the carbon-rich 
late-type star IRC+10216, where their similar distributions are a puzzle for 



photochemical models (Il6l ). In cometary atmospheres, species as complex as 



C4H have been inferred, whose presence is not readily explained by the stan- 



dard parent-daughter photochemical models (Il7l ). 



Many theoretical studies of the electronic structure of small carbohydride 



molecules exist in the literature (1181; Il9l; l20t l2ll; |22|; l23l. e.g.,). However, most of 



them are limited to the ground and lowest few excited electronic states. Only 
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few studies of higher states exist, most of them performe d by Peyerimhoff and 



her co-workers starting more than 30 years ago, e.g., (1241 ; l25f ) As will be shown 
here, those results are still highly relevant and provide an excellent starting 
point for further studies. 

The aim of this paper is to provide insight into the photodissociation processes 
of carbohydride molecules through ah initio quantum chemical calculations of 
the vertical excitation energies and oscillator strengths for as large a number 
of excited electronic states as feasible with current programs. By assuming 
that most of the absorptions with energies above the dissociation limit lead 
to destruction, estimates of the (upper limits of) photodissociation rates un- 
der interstellar conditions can be obtained. Even though such calculations are 
necessarily limited to the smaller members of the carbo-hydride family, they 
do provide quantitative constraints to test the basic interstellar chemical net- 
works. Specifically, C 2 H, C 3 , C 4 , I— and c— C 3 H, I— and c— C 3 H 2 , HC 3 H, 
I— C 4 H and I— C 5 H are studied here. Except for HC 3 H, all these species have 
been chosen to have been observed in interstellar space. 



2 Methods 



2.1 Electronic structure calculations 



For small diatomic and triatomic molecules, quantum chemical calculations of 
the potential energy surfaces and transition dipole moments combined with 
dynamical calculations of the nuclear motions can provide photodissociation 
cross sections and oscillator strengths that agree with experiments to better 



than 20-30%, e.g., (J26J; |27|; [28|) (see (I29l ) for an early review). For the poly- 



atomic molecules considered here, calculations of the full potential surfaces 
including all degrees of freedom become very time consuming, as do the multi- 
dimensional dynamics. Moreover, such detail is not needed in order to compute 
photorates, since those are largely determined by the potentials and transition 
moments in the Franck-Condon region. The simplest alternative is therefore 
to only compute the vertical excitation energies and transition dipole moment 
of the molecule at its equilibrium position and assume a certain probability 
that absorption into each excited state above the dissociation limit leads to 
dissociation. For the examples cited above, such an approach leads to similar 
rates within the accuracy of the calculations. The focus of our calculations is 
therefore on electric dipole-allowed transitions to states lying above the low- 
est dissociation limit but below 13.6 eV. The interstellar radiation field has 
a broad spectrum from visible to extreme ultraviolet wavelengths, peaking in 
intensity around 7 eV and cutting off at the atomic H ionization potential at 
13.6 eV (Br 
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All calculations presented here were performed with the Wuppertal-Bonn 
MRDCI set of p rog rams as implemented in the GAMES S-UK program pack- 
age version 7.0 (13 ll ). For C, the TZVP atomic orbital basis set was used, and 
for H, the DZP basis set ((32). To allow for a proper description of molecular 
Rydberg states, two diffuse p and two diffuse d functions were put on specific 
sites. For C2 and C2H, this was the middle of the CC bond; for C2 and C4H 
the middle of the C1-C2 and C3-C4 bonds; for C 3 , /-C 3 H, /-C 3 H 2 and HC 3 H 
the middle C; for C 5 on C2 and C4; for c— C 3 H on the C connected to H; and 
for c— C 3 H 2 on the lone C. These sites were chosen to avoid linear dependen- 
cies in the atomic orbital basis set due to large overlap of diffuse functions 
positioned on adjacent atoms. Cartesian d functions were used, of which the 
spherical component served as the 3s basis function. 



Molecular orbitals were generated using up to the maximum number of ref- 
erence states, 255. The selection threshold was generally set at 0.5 /xHartree. 
The total number of configurations included in the configuration interaction 
(CI) calculation ranged from 200000-400000 per symmetry. The sum over 
the coefficients-squared in the final CI wave function is typically 0.95 for the 
smaller species, dropping to 0.9 for the larger molecules. The CI energies were 
extrapolated and corrected using the Davidson extrapolation (1331 ). The aim 
was to compute as many excited electronic states as feasible, up to 9 per sym- 
metry. For the lower states, comparison with existing calculations and experi- 
ments indicates accuracies within 0.3 eV, generally better. Oscillator strengths 
to the lower states agree within 30% or better. For the higher states, typically 
the 5'th root and higher per symmetry, the accuracy decreases because many 
states and orbitals can mix. Nevertheless, such calculations should still pro- 
vide insight into the location of those states, in particular whether they are 
above or below the ionization potential and below the 13.6 eV cutoff of the 
interstellar radiation field. Moreover, the magnitude of the oscillator strengths 
(strong or weak) should be reliable. Note that the precise values of the excita- 
tion energies are not so important for the purposes of calculating interstellar 
photorates because of the broad range of incident energies. The only exception 
is possible overlap with Lyman a radiation at 10.2 eV, which is important for 
certain astrophysical environments (1341 ) . 



In all calculations, the largest Abelian subgroup of the full group, C 2v , 
was used, with the molecule put along the z— axis. The doubly degenerate 
states are then split into the B\ and B 2 irreducible representations for II and 
$ states, and into A\ and A 2 for the A states. In general, the degeneracies for 
the A states are recovered within a few hundreds of an eV in the calculations. 
The average of the two A\ and A 2 values is tabulated here. The A\ irreducible 
representation also contains the S + states, and the A 2 the S~ states. Some 
excited states with electronic angular momentum higher than 2 (in particular 
$ states) are also found in our calculations but they are not tabulated here 
since the electric dipole transition moments from the ground state to these 
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states are zero, implying that they do not contribute to photodissociation. The 
valence or Rydberg character of the states is determined from calculations of 
the x 2 , y 2 and z 2 expectation values. 



Calculations have been limited to the equilibrium geometry corresponding to 
the lowest energy. For C 3 , C 2 H, C 4 , C 4 H and C 5 H this is the linear geometry. 
C3H and C3H2 have been detected in various isomeric forms so both cyclic 
and (near-)linear forms have been studied. 

The dissociation energies for C n H are calculated as the difference between 
the ground state energy of C n H at its equilibrium geometry and the energy of 
ground state C n , also at its equilibrium geometry, calculated in the presence of 
an H atom positioned at a distance of 20 Bohr from the center of mass. In this 
way the errors associated with the lack of size consistency are significantly 
reduced. Note that the spatial symmetry of the molecule and the fragment 
generally differs. 



2.2 Photodissociation rates 

Oscillator strengths to all excited states have been calculated according to 



where all quantities are in atomic units (a.u.), fi u i is the transition dipole 
moment from lower state I to upper state it, and g u i is a degeneracy factor 
which is 2 for a II <— £ transitions and 1 for any other transition. For linear 
molecules, the dipole moment operator contained in \i u i corresponds to — X Zj 
for transitions between states with the same value of the electronic angular 
momentum projection quantum number A and to — XX x j + Wj) I for tran- 
sitions with A A = ±1. The computed excitation energies AE u i were used in 
this formula, not any experimental values. Hence, differences with other work 
can stem from differences in both fi u i and in AE u i. 

The photodissociation rate of a molecule can be computed from 



where a is the photodissociation cross section in cm 2 and I is the mean 
intensity of the radiation in photons cm -2 s _1 A -1 (34). Under interstellar 
conditions, only single photon processes are important. For photodissociation 




- g ui AE ul fi ul 



(1) 




(2) 
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initiated by line absorptions (e.g., predissociation), the rate becomes 



ire 



2 



k pT = —^KifuiriJui s- 1 (3) 
1 mc z 



where r] u the dissociation efficiency of state u, which lies between and 1. The 
numerical value of the factor Tce 2 /mc 2 is 8.85 x 10 -21 in the adopted units with 
A in A. The total photodissociation rate of a molecule is obtained by summing 
over all channels. In this work, no dynamical calculations are performed to 
obtain continuous cross sections for dissociative states. Hence, all photorates 
are computed with Eq. (3). It is furthermore assumed that the dissociation 
efficiency T] u —1, either through direct dissociation in a repulsive state or by 
predissociation. The motivation for this choice is that for larger molecules, 
internal conversion to a lower (dissociative) electronic state is usually much 
more rapid than any radiative decay rates owing to the high density of states 
(jiil ). Specific experimental evidence of high dissociation efficiencies will be 
presented in the results section for individual molecules. 



3 Results 



3.1 C 2 H 



The C2H radical, detected in interstellar clouds since 1974 (Q]), is an impor- 
tant step in the formation of longer carbon chains. In comets, it could be 
a photodissociation product of C2H2 and a precursor of the widely observed 
C2 molecule. Given also its importance in combustion processes, this radical 
has received ample theoretical attention, starting with the papers by Shih, 
Peyerimhoff and co-workers (1241 ; l36l ) and culminating with the more recent 
two-dimensional potential surfaces by Duflot et al. (1371 ). 



As a test of our computational procedure, we present in Table 1 our computed 
vertical excitation energies, oscillator strengths and the main configurations 
of each state. Our energies generally agree well within 0.3 eV with those of 
Koures and Harding (1381 ) (their CI + DV2 results) and lie in between those 
of Duflot et al. (1371 ) and Shih et al. (24). Close inspection shows that the 
larger differences are usually for states with Rydberg character (which are 
not included in Duflot et al. (1371)) or for states with a strong interaction with 
a neighboring state, leading sometimes even to a switch in the character of 
the state. For example, the 4 2 £ + state in our calculation has Rydberg 3p n 
character whereas it has ... 4cr 1 l7r 3 5<T 2 37r 1 in Duflot et al., resulting also in a 
very different transition dipole moment (see below). 
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Our computed A 2 U — X 2 X + transition dipole moment of 0.22 a.u. is very 
close to that of 0.23 a.u. computed by Duflot et al. (37) and Peric et al. (j39l ). 
Because of the slightly higher excitation energy in our work, the oscillator 
strengths show somewhat larger differences. Comparison for other states is 
difficult since Duflot et al. do not give any numerical values. However, their 
Figure 7 shows that neither the 3 and 4 2 £ + , nor the 2 and 3 2 IT states 
have significant transition dipole moment. This is generally consistent with 
our results. For the 5 2 £ + state we find a huge transition dipole moment of 
1.1 a.u. due to the Rydberg character of our wavefunction. Thus, this state 
around 10 eV will dominate the interstellar C2H photodissociation but the 
higher 2 n states in the 8.5-10.5 eV range can also contribute significantly up 
to the ionization potential of ~11.4 eV (Table 11). 



3.2 C 3 



The C3 molecule was detected in cometary spectra in 1882 through its A — X 
system at 4050 A (J40j) and in interstellar clouds more than a century later 
(41; [sl). The molecule is also seen in the atmospheres of cool carbon stars 



through its mid-infrared (1421 ) and far-infrared (1431 ) transitions (see (1231 ) for 
review) . 

C3 is a linear molecule with a ground X 1 S+ state. Its dissociation energy is 
computed to be about 4.6 eV, whereas its ionization energy is around 12 eV 
(Table 11). Quantum chemical studies range from the early work by Cha- 



balowski et al. (144l ) to the recent calculations by Terentyev et al. (145l ). Table 2 
summarizes our computed excitation energies, together with the oscillator 
strengths and main configurations. For the low-lying valence states, our re- 
sult s ag ree to better than 0.3 eV with the MR-AQCC values of Monninger et 
al. (1461 ) . These results also show that the X E+ <— X l Y,+ (\n u — > ln g ) tran- 



sition at 8.17 eV has by far the largest oscillator strength, as predicted first 
by Pitzer and Clementi (1471 ). This is confirmed by the combined experimental 



and theoretical study by Monninger et al. (46j) whose 1100-5600 A spectrum 
demonstrates that the <— X l T^ transition around 1700 A is indeed the 
strongest band. Their experiments in Ne and Ar matrices show a broad band 
but with some progressions superposed, which are evidence for interstate vi- 
bronic coupling with adjacent H g states. Indeed, the state is predicted 
to be unstable to bending, leading to avoided crossings with the lower-lying 
Tig states (|46l ). These interactions can also lead to dissociative channels to 
produce C2 + C. 
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3.3 C 3 H 



Both linear and cyclic C3H have been detected in the interstellar medium by 
their transitions at millimeter wavelengths (|48l ; l49l ). The 2 U ground state of 
I— C 3 H lies about 0.4-0.6 eV above the 2 B 2 ground state of c— C 3 H (j5ol ). with 
a small barrier toward the cyclic state. Hence, when C3H is produced in its 
linear form by some sequence of ion molecule or other chemical processes, it 
can be stable, and both isomers are therefore considered in this work. 



3.3.1 I-C3H 

The geometry of the 1 2 I1 state was taken from the experimental work of Mc- 
Carthy and Thaddeus ([69). Table 3 summarizes our computed vertical ex- 
citation energies, oscillator strengths, and the corresponding configurations. 
Comparison with the CASSCF results of Ding et al. (51) shows good agree- 
ment for the lowest states. The dissociation energy of I— C3H to C3 + H is 
computed at 3.3 eV. Of the dip ole- allowed states above this dissociation limit, 
the higher 2 I1 states around 7.8 eV have the largest oscillator strengths, in 
particular the 8 2 I1 state, which will dominate the photodissociation of the 
molecule. 



3.3.2 c-C 3 H 

The c— C3H radical has C 2v geometry in its 1 2 B 2 ground state, with the 
equilibrium structure for our calculations taken from Yamamoto and Saito 



(I52l ) . Table 4 summarizes our computed vertical excitation energies, oscillator 



strengths and configurations. Comparison with the CASSCF results of Ding 



et al. (15 ll ) shows again good agreement for the lowest states to 0.1-0.2 eV. 
The dissociation energy of c— C3H to C3 + H is around 4.3 eV whereas its 
ionization potential is computed to lie at 9.6 eV (Table 11). Of the dipole- 
allowed states above the dissociation limit, the 2 2 Ai state around 5 eV and 
the higher 2 A 2 states around 7.5 eV have the largest oscillator strengths, but 
the sum over the other states is comparable. Thus, there appear to be many 
potential routes to photodissociation for c— C3H. 



3.4 C 3 H 2 



Cyclopropenylidene, c— C3H2, was the first cyclic molecule to be detected in 
the interstellar medium (jj) and subsequently found to be ubiquitous through- 
out the Galaxy, even in diffuse gas (J53J; |9|). The near-linear form HCCCH 
(propargylene) , denoted here as HC 3 H, lies 0.4-0.8 eV higher in energy (154 ; 
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50|), depending whether the zero-point vibrational energy is included. Since 
it has a near-zero dipole moment, it has not yet been observed in interstellar 
clouds through radio transitions but is likely present as well. There are various 
other stable isomers of C3H2, of which another linear form, H 2 CCC (vinyli- 
denecarbene or propadienylidene) , denoted here as I— C3H2, was discovered in 
interstellar space in 1991 (I55l ). It lies ~0.6 eV above the cyclic ground state in- 



cluding zero-point vibrational energy corrections (154 ) . Hence, all three isomers 



are considered in this work. Although there are many theoretical studies of the 



ground states of the various C 3 H 2 isomers dating back to 1976 (1561 ; Il8t 1571 ). 
the excited states are largely unexplored. 

The equilibrium structures of the electronic ground states are taken from 
Seburg et al., their Figure 5 (154 ). Our calculated relative ordering of the C3H2 



isomers is consistent with previous findings but due to different amounts of 
recovered correlation energy, our energy splittings are larger than in other 
studies. 

The adiabatic ionization potentials of the various C3H2 isomers lie at 8.96 
(HC 3 H), 9.2 (c-C 3 H 2 ) and 10.4 eV (Z-C 3 H 2 ), respectively (Table 11). Above 
these thresholds, the photoionization efficiency of all three C3H2 isomers in- 
creases rapidly so that photoionization will become the main pathway (jslj ). 
Some of these ionizations are likely to be dissociative, but this option is not 
considered here. 



3.4.I HC 3 H 

The HC3H radical has C 2 geometry with a 3 B ground state. In this case, 
the C2 axis was put through the middle C atom and the middle of the 
H1-H2 line. Table 5 summarizes our computed vertical excitation energies, 
oscillator strengths and configurations. For the lowest four excited states, 
our results are in good agreement with those of Mebel et al. (I59I ) (their 
MRCI+D(4,8)/ANO(2+) results). The dissociation energy of HC 3 H to Z-C 3 H 
+ H is computed to be around 3.1 eV, with the dissociation energy to C 3 + 



H2 perhaps even less (1591). Of the dipole- allowed states above the dissociation 
and below the ionization limit at 8.96 eV, the higher 3 A states at 6.2 and 7.5 
eV, both of which have n — > 7r* character, have the largest oscillator strengths. 



3.4.2 c-C 3 H 2 

The c— C 3 H 2 molecule has a 1 A 1 ground state in its C 2v geometry, with a large 
dipole moment of 3.4 Debye owing to the two unpaired electrons on one of 
the three carbon atoms. Table 6 summarizes our computed vertical excitation 
energies and oscillator strengths, which again agree well in terms of energies 
with those of Mebel et al. (159l ) . although our oscillator strengths are some- 
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what lower. Compared with c— C3H, c— C3H2 has only few low-lying electronic 
states, and all of the dip ole- allowed transitions lie above the dissociation en- 
ergy of 4.4 eV (59|). Those to the Q 1 A 1 (tt tt*) at 9.4, 7 1 B 1 at 10.6 and 
9 l Bi state at 11.2 eV have the largest oscillator strengths. However, none of 
these states lie below the ionization potential of c— C3H2 at 9.2 eV. Thus, its 
photodissociation rate will differ substantially whether or not these states are 
included (see §4). 



3.4.3 l-C 3 H 2 



The vinylidenecarbene or propadienylidene isomer of C3H2, denoted here as 
/— C 3 H 2 , also has 1 A 1 symmetry in its C 2v ground state geometry. Our com- 
puted vertical excitation energies (Table 7) of 1.86, 2.44 and 5.54 eV to the 
A 1 A 2 , B 1 Bi and C(2) l Ai states are consistent with the experimental adiabatic 
values (J6Cj) of 1.73, 2.00 and 4.84 eV, respectively, and agree within 0.1-0.2 eV 
with the vertical values computed by Mebel et al. (J59j) The dissociation energy 
of I— C 3 H 2 to C 3 H + H is computed to lie around 3.9 eV, so that only the 
C(2) 1 Ai and higher states can lead to photodissociation. The C state shows 
a well-resolved progression in its electronic spectrum (60), but the resolution 
of those data is not high enough to measure predissociation rates. Thus, we 
consider the photodissociation both with and without taking the C state into 
account. 



Of the dipole- allowed transitions below the ionization potential at 10.4 eV, 
those to the higher 1 Ai states around 9 eV have the strongest oscillator 
strengths. These states all have 2b 2 — > n*b 2 character and likely belong to 
the Rydberg series converging to the lowest ionization potential. Similarly, 
the higher l B 2 states belong to the 2b 2 — > nai Rydberg series. 



3.5 C 4 



The lowest energy ^ 3 S fl ground state of C 4 occurs for the linear geometry. 
C 4 has been searched for in diffuse interstellar clouds through its 3 £~ ^ 3 S g 
transition around 3789 A but not yet detected (16 ll ). However, a pattern of 
bands at 57 um observed with the Infrared Space Observatory toward a handful 
of objects (j62j) is consistent with transitions in the z/ 5 bending mode (1631 ) . 
Rhombic C 4 , which is almost isoenergetic with I—C4, has not yet been detected 
in interstellar space. 



Our calculations use the I— C 4 equilibrium values computed by Botschwina 
(J64j), which are consistent with experiments (1231) . The results are presented 



in Table 8. Our excitation energies of the lower states are in goo d agreement 
with experiments and with the calculations of Masso et al. (65J), except for 
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the 2 3 A U state. Note that the experimental values refer to To rather than T e , 
which can differ by a few tenths of eV. Our computed dissociation energy D e 
is 4.7 eV, in good agreement with the experimentally inferred D value of 4.71 



±0.15 eV (J66|). 



The lowest excited state above the dissociation limit to which dipole allowed 
transitions are possible is the 3 3 U U state at 5.1 eV. Indeed, photofragment yield 
spectra in the 2.22-5.40 eV range by Choi et al. (66) show significant single- 
photon dissociation into C3 + C at >5.2 eV, with a minor channel to C2 + 
C2. Their spectra are well reproduced by phase space theory models in which 
the product state distributions are statistical. This implies that absorption 
into the excited state is most likely followed by rapid internal conversion to 
the ground state potential energy surface with no barriers present along the 
dissociation coordinates. 

By far the strongest absorption occurs into the 2 3 S~ state around 6.95 eV, 
which has an oscillator strength of 1.56. Even though higher S~ and U u states 
will contribute, the 2 3 S~ channel will dominate the interstellar photodissoci- 
ation of C4, if indeed every absorption is followed by dissociation. 



3.6 l-C 4 H 



The I— C4H (butadiynyl) radical was detected in the envelopes of carbon-rich 



evolved stars nearly 30 years ago (1671 ). It was subsequently found to be very 



abundant in cold dark clouds like TMC-1 (1681 ) and even detected in comets 



(17). It is one of the carbon-bearing molecules found at the edges of PDRs 
(115l). where it can photodissociate into smaller species. 



C4H has linear symmetry with a S + ground state. The equilibrium coordi- 



nates in our calculations were taken from Ref. (1691 ) . Table 9 summarizes our 
computed vertical excitation energies, oscillator strengths and configurations. 
Comparison with the results of Graf et al. 0) shows good agreement in both 
energies and transition dipole moments for the lowest states, but poor agree- 
ment for the higher 2 IT states, where the Graf et al. energies are generally 
lower by up to 1 eV. There are two main resons for this. First, Graf et al. 
did not include diffuse (Rydberg) functions in their basis set, which start to 
become important for the higher states. Second, their CASPT2 perturbation 
method does not guarantee a lower bound to the energies. Our calculations 
agree in the fact that none of the higher lying 2 II states have large oscillator 
strengths. 

The threshold for photodissociation, corresponding to dissociation of C4H to 
C4 + H, lies around 4.7 eV (Table 11). Of the dipole-allowed transitions below 
the ionization potential at 9.6 eV, the 4 and 5 2 S + states at 7.7 and 8.7 
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eV, respectively, have orders of magnitude larger oscillator strengths than 
other states and will thus dominate the photodissociation. Both states have 
l7r — > rnr Rydberg character. Thus, our overall photodissociation rate will be 
significantly larger than that using the data from Graf et al. (1701 ) 



3.7 l-C 5 H 



Like C 4 H, /— C 5 H is detected toward the carbon-rich evolved star IRC+10216 
(fnl ) and in cold dark clouds (72). It also has linear symmetry but with a 
2 I1 ground state. The equilibrium coordinates in our calculations were taken 



from Ref. (]69l ). Table 10 summarizes our results. Comparison with Haubrich 
et al. (|73l ) shows excellent agreement in both energies and oscillator strengths, 
except for the higher-lying 2 E + states. This difference is likely due to the ex- 
plicit inclusion of Rydberg states in our work. For the 3 and 4 2 U states, the 
oscillator strengths differ by a factor of 2 but these states show considerable 
interaction between the 2n — > 37r and 3n — > An (n — ► 7r*) excitations. Differ- 
ences in mixing ratios can lead to large changes in transition dipole moments 
to individual states, but not in energies. The computed excitation energies 
for the lowest two excited states also agree well with those measured and 
computed by Ding et al. (1741 ) . 

The threshold for photodissociation lies around 3.6 eV (Table 11). Of the 
dipole-allowed transitions below the ionization potential at 7.4 eV, both the 4 
2 n (2tt -> 3tt) state at 4.2 eV and the 6 2 IT (1 7r — > 3tc) state around 6.1 eV have 
large oscillator strengths. All other states have typical oscillator strengths of 
a few xlO -3 and thus contribute at a lower level. 



4 Interstellar photodissociation rates 



In Table 12, the computed photodissociation rates in the unshielded interstel- 
lar radiation field cf. Draine (j3p| ) are presented, using the oscillator strengths 
given in Tables 1 to 10 and asuming i] u =l for all states. Thus, these rates 
should be regarded as upper limits. Only states above the dissociation limit 
and below the ionization potential have been taken into account (see Table 11 
for adopted values). No corrections have been made for possible higher-lying 
states below the ionization limit not computed in this work. Since the oscilla- 
tor strengths for Rydberg states decrease roughly as 1/n 3 , it is assumed that 
any such corrections would be small since the lowest Rydberg members are 
calculated explicitly. For reference, inclusion of a hypothetical state at 9 eV 
with an oscillator strength of 0.1 would increase the photodissociation rates 
by only 3.5 x 10" 10 s _1 . 
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For C2H, our new rate is a factor of 3 larger than that given in van Dishoeck et 



al. (134j ). which was based on the energies and oscillator strengths of Shih et al. 
(1241 ). The increase is mostly due to the higher 2 E + and 2 I1 states which were 
not computed in that work. For C3, the new rate is only 30% larger, mostly 
because the photodissociation rate of this molecule is dominated by the very 
strong absorption into the l 1 ^ state around 8 eV, which was included in 
previous estimates. 

It is seen that the rates for the various carbon-bearing molecules span a range 
of a factor of 8, with the rates for the bare carbon chains (C3, C4) being largest 
and those for the odd- numbered C n H species lowest. Of the different C 3 H 2 iso- 
mers, I— C3H2 has the largest photodissociation rate and c— C3H2 the smallest 
by a factor of 3. However, c— C3H2 differs from the other isomers in that it 
has several states with large oscillator strengths above the ionization potential 
of 9.15 eV. If those states were included, the c— C3H2 photodissociation rate 
would be increased by a factor of 2. For I— C3H2, the rate would drop from 
5.1 x 10~ 9 to 4.1 x 10~ 9 s" 1 if the C state is not included. 



In spite of the range of values, all rates are above 10~ 9 s , coresponding to a 
lifetime of less than 30 yr at the edge of an interstellar cloud. In regions such 
as the Orion Bar and the Horsehead nebula, the radiation field is enhanced by 
factors of 10 3 — 10 5 compared with the standard field adopted here, decreasing 
the lifetimes to less than 1 month. Thus, there must be rapid production routes 
of these molecules in order to explain their high abundances in UV-exposed 
regions. This conclusion is not changed if only 10%, say, of the absorptions 
would lead to dissociation rather than the 100% assumed here. As argued in 
§2 and 3, it is plausible that a substantial fraction of the absorptions lead to 
dissociation for these larger molecules so that the upper limits should be close 
to the actual values. 



5 Conclusions 



We have presented vertical excitation energies, transition dipole moments and 
oscillator strengths for states up to the 9'th root of each symmetry for several 
carbonaceous molecules of astrophysical interest. For lower- lying states, good 
agreement is generally found with previous studies. Several new, higher-lying 
states with large oscillator strengths, often of Rydberg character, are revealed 
in this work. 

The calculated photodissociation rates of the small carbon-bearing molecules 
studied here are substantial, leading to lifetimes at the edges of interstellar 
clouds of less than 30 yr. These high rates assume that all absorptions above 
the dissociation limit indeed lead to dissociation, so that the rates should be 
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viewed as upper limits. Further experimental work is needed to quantify the 
dissociation efficiencies for the strongest states found in this work. 
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Table 1 

Vertical excitation energies, oscillator strengths and dominant configurations for 
C2H at the ground state equilibrium geometry 







Energy (eV) 




f el 


Dominant 


Type" 


State 


This 


Ref. 


Ref. 


Ref. 


This 


Configuration 






work 


(38) 


(2D 


(M) 


work 






l 2 £+ b 


0.00 


0.00 


0.00 


0.00 




..Aa^aHir 4 


V 


2 2 S+ 


7.06 


6.73 


6.63 


7.32 


4.0(-4) c 


...4cj 2 5cJ 1 l7r 3 27r 1 


V 


3 2 £+ 


8.63 


8.11 


8.19 


9.60 


4.0(-5) 


..Aa^aHir 4 


V 


4 2 S+ 


9.28 


9.01 


9.21 


9.18 


3.0(-3) 


...4o- 2 5a 1 l7T 3 3px 1 


R 


5 2 £+ 


10.09 


10.09 






2.8(-l) 


..Acr 2 5a 1 lir 3 4px 1 


R 


6 2 £+ 


10.28 








1.0(-5) 


..Aa 2 ha x l%^px l 


R 


i 2 n 


0.68 


0.60 


0.54 


0.96 


1.7(-3) 


...4fj 2 5o- 2 lvr 3 


V 


2 2 n 


7.63 


7.29 


7.07 


8.11 


1.0(-2) 


...4cj 2 5cJ 2 l7r 2 27r 1 


V 


3 2 n 


8.39 


8.17 


8.05 


9.96 


1.0(-1) 




R 


4 2 n 


9.00 


8.68 


8.34 


8.48 


3.0(-2) 


...4o- 2 5a 1 l7T 3 3pz 1 


R 


5 2 n 


9.47 




8.70 




3.0(-2) 


...4cr 2 5(T 1 l7r 3 3s 1 


R 


6 2 n 


9.96 




8.80 


9.22 


2.2(-2) 


...4cr 2 5cj 1 l7r 3 4^ 1 
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7 2 n 


10.06 




9.25 
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l.l(-2) 
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8 2 n 
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R 
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8.25 


8.12 
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V=Valence; R=Rydberg; M=Mixed in this and subsequent tables 



b Ground state energy including Davidson correction: -76.432130 Hartree 
c Notation x(-y) in this and subsequent tables indicates x x 10~ y 
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Table 2 

Vertical excitation energies, oscillator strengths and dominant configurations for C3 
at the ground state equilibrium geometry 





Energy (eV) 
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Ref. 
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a Ground state energy including Davidson correction: -113.735304 Hartree 
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Table 3 

Vertical excitation energies, oscillator strengths and dominant configurations for 
I— C3H at the ground state equilibrium geometry 





Energy (eV) 
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a Ground state energy including Davidson correction: -114.381395 Hartree 
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Table 4 

Vertical excitation energies, oscillator strengths and dominant configurations for 
c— C3H at the ground state equilibrium geometry 





Energy (eV) 
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1.8(-2) 


Qcl\ — ► 




V 


3 2 A 2 


C CO 




4.2(-4) 


6czi — > 




v 


4 2 A 2 


7 




2.7(-2) 


(6a!) 2 


— * 3b 2 2a 2 


R 


5 2 A 2 


7 ^4 

i .04: 




1.8(-2) 


6ai — > 


Abi 


R 


6 2 A 2 


O.fJO 




1.3(-2) 


5oi — > 


2bi 


R 


7 2 Ao 


8.72 




9.6(-3) 


36 2 -» 


la 2 


R 


8 2 A 2 


8.99 




7.3(-3) 


6ai — > 


Sbi 


R 


9 2 A 2 


9.46 




4.0(-3) 


6ai — > 


46i 


R 


l 2 B l 


3.47 


3.49 




3b 2 -» 


26i 


V 


2 2 B l 


4.70 


4.63 


21 


36 2 -> 


36i 


V 


3 2 B 1 


5.73 






601 — > 


la 2 


V 


a2 td 
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Table 5 

Vertical excitation energies, oscillator strengths and dominant configurations for 
HC3H at the ground state equilibrium geometry 

Energy (eV) f el Dominant Type 
State This This Configuration 

work work 



1 3 ™ 

1 D 


n nn 
u.uu 






...DO Id 60 40 


V 


93 R 
Z D 


A 17 
4.0/ 


1 8/ 


/i 


oO — ' 40 


V 


^3 R 


A ^8 


1 01 


-0 j 


DO — ' 1 a 


V 


/l3 R 
4 D 


O.D1 


/I 9/ 
4.ZI 


-0 ) 


7„ . Q„ 

1 a — ► ya 


rt 


^3 R 


o.oy 


1 .01 




7„ . Q„ 

/ ci — > oa 


v 

V 


6 3 B 


6.08 


6.7( 


-4) 


46 -> 56 


V 


T A B 


6.67 


7.0( 


-6) 


7a -► 10a 


R 


8 3 B 


6.93 


6.8( 


-4) 


46^66 


R 


9 3 B 


6.97 


1.9( 


-4) 


7a -> 11a 


R 


1 3 A 


4.23 


1.7( 


-2) 


7a -»• 46 


V 


2 3 A 


4.33 


2.4( 


-4) 


36 -► 7a 


V 


3 3 A 


5.66 


9.3( 


-4) 


46 -► 8a 


V 


4?A 


5.97 


1.7( 


-3) 


7a -► 56 


V 


5 3 A 


6.29 


1.5( 


-1) 


7a -► 66 


R 


6 3 A 


6.68 


1.2( 


-3) 


46^ 9a 


R 


7 3 A 


6.95 


4.8( 


-3) 


46 -► 10a 


R 


8 3 A 


7.52 


2.4( 


-1) 


46 -► 12a 


R 


9 3 A 


7.73 


2.2( 


-3) 


46 -► 13a 


R 



a Ground state energy including Davidson correction: -114.924751 Hartree 
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Table 6 

Vertical excitation energies, oscillator strengths and dominant configurations for 
c— C3H2 at the ground state equilibrium geometry 



State 



Energy (eV) 
This 



f el 

This 



Dominant 
Configuration 



Type 





work 


work 








0.00 






...5af6af 16f 26|36| 


V 


2 l Ai 


6.31 


1.1 


^7 


6ai — > 8ai 


R 


3 1 A 1 


6.85 


6 71 


L ) 


6ai — ► 9ai 


R 


A 1 A 1 


7.82 


1.1 


-2\ 


6ai — > 10ai 


R 


5 1 A 1 


8.11 


2.9( 


-3) 

°7 


601 — > 9ai 


R 


„ 1 . 

6 1 A 1 


8.57 


7.0( 


-2) 
^7 


6ai — > llai 


R 


7 1 A 1 


9.37 


1.2 


-1) 


lfei -> 26i 


V 


8 1 A 1 


9.70 


4 8( 


-2) 


6ai — > 12ai 


R 


_ 1 . 


10.94 


8.7( 


-4) 

^7 


6ai — > 12ai 


R 




4.89 


2.8( 


-2) 


6ai — > 56i 


V 


2 1 .Bi 


6.65 


2.4( 


-3) 

J 7 


36 2 — ► la 2 


V 


3 1 B 1 


6.91 


1.9( 


-2) 
^7 


6ai — > 2b\ 


R 


. 1 T 

4 1 B 1 


8.25 


6.3( 


-3) 


6ai — > 36i 


R 


_ 1 _ 


8.63 


1.1 


-2) 


Ibi -»■ 8ai 


R 




10.30 


4.3( 


-2) 

^7 


lfei -> 9ai 


R 


7 1 J B i 


11.16 


2.6' 


-2\ 


36 2 — ► 1«2 


R 




11.53 


1.8' 


-4) 

^7 


I61 -> 10ai 


R 


9^1 


12.30 


1.3( 


-3) 


601 -► 46i 


R 




6.68 


1.4i 


-4) 

^7 


6ai — ► 662 


R 


2 1 5 2 


8.13 


4.4( 


-3) 

°7 


601 -► 76 2 


R 




8.29 


1.1( 


-2) 


6a 1 — > 4&2 


V 


I 1 „ 

4 1 £ 2 


9.01 


6 6( 


-2) 


l&i — ► la 2 


R 


5^2 


9.20 


5.8( 


-3) 


362 — ► 8ai 


R 


6^2 


9.99 


3.0( 


-2) 


362 — ► 9ai 


R 


7^2 


10.58 


3.0( 


-1) 


6ai — > 5&2 


V 


8^2 


11.03 


1.6( 


-5) 


3&2 - ► 8ai 


R 


9^2 


11.23 


1.0( 


-1) 


6ai — > 862 


R 


1% 


3.92 






6ai — > la2 


V 


2% 


7.03 




23 


3&2 — ► 56i 


V 


3% 


8.46 






6ai — > 202 


R 


/ll /t_ 


n ok 











Table 7 

Vertical excitation energies, oscillator strengths and dominant configurations for 
/ — C3H2 at the ground state equilibrium geometry 



State 



Energy (eV) f el 
This This 



Dominant 
Configuration 



Type 





work 


work 








0.00 








V 


2 1 A 1 


5.54 


1 3 


1 \ 

/ 


I61 26i 


V 




5.98 


4.7 


-2) 


(26 2 ) 2 - (2&!) 2 


V 


A 1 A 


7.79 


9 8 


_2) 


26 2 36 2 


R 


5^i 


8.87 


1.5 


-2) 


26 2 -» 46 2 


R 


^ A 


9.00 


3 


-1) 


26 2 56 2 


R 


7 1 Ai 


9.38 


3.0 


-1) 


26 2 -»• 662 


R 


— 1 A 

8^i 


9.84 




_2) 


26 2 76 2 


R 


9^i 


10.02 




_2) 


I61 -► 36i 


R 




2.44 


O.O 


-3") 


7ai -> 26i 


V 


2 1 #i 


6.36 


1.3 


-2) 


7oil6i -► (26i) 2 


V 


3 1 Bi 


8.01 


1.7 


_2) 


7a 1 — > 36i 


R 


4 1 #i 


8.82 


1.3 


_3) 


7ai26 2 -> 26i56 2 


V 




8.86 


3.8 


.3) 


7ai -> 46i 


R 




8.88 


7.9 


_2) 


\b x -> 10ai 


R 


7 1 J B i 


9.11 


2.1 


_2) 


7a 1 — > 56i 


R 




9.85 


3.6 


_3) 


7a 1 — ► 661 


R 




10.13 


1.9 


_2) 


I61 -► 8ai 


R 


l l B 2 


6.06 


2.4 


-4) 


6ai26 2 -► (26i) 2 


V 


2 1 B 2 


6.92 


1 3 


/ 


26 2 -> 8ai 


V 


3 X B 2 


8.07 


4.0 


-2) 


26 2 -> 9ai 


R 


4}B 2 


8.87 


4.2 


-4) 


26 2 -> 10ai 


R 


b 1 B 2 


8.98 


2.0 


-4) 


26 2 -► Hai 


R 


&B 2 


9.36 


1.9 


-2) 


26 2 12ai 


R 


1 X B 2 


9.68 


7.2 


-5) 


26 2 -► 13ai 


R 


8 X B 2 


10.14 


3.4 


-2) 


26 2 -> 14ai 


R 


9 X B 2 


10.50 


<1 


'-6) 


26 2 -> 15ai 


R 


l l A 2 


1.86 






26 2 -► 26i 


V 


2 l A 2 


6.68 


24 


16i26 2 (260 2 


V 


&A 2 


7.71 






26 2 36i 


R 


Al a. 


Qn 






1f,_ V /If,. 


\r 



Table 8 

Vertical excitation energies and oscillator strengths and dominant configurations for 
I— C4 at the ground state equilibrium geometry 





Energy (eV) 


f el 


Dominant 


Type 


State 


This 


Ref. 


Exp. 


This 


Configuration 






work 


(65) 


(65) 


work 






l 3 £" a 

y 


0.00 


0.00 


0.00 




..Aal5* 2 g lTr£lTr 2 g 

u, y u, y 


V 


2 3 S" 
9 


6.40 


6.11 






...4^5cj^l7r 3 l7r 3 

u g u g 


V 


3 3 £" 
9 


6.79 








ll y ll y ll 


V 


4 3 S" 
9 


7.05 








..A0 2 Ji0 2 \Ttl\K 2 2'K\ 

u y u y u 


V 


1 3 S+ 
9 


5.80 








LL y LL y LL 


V 




3.61 


3.74 


3.27 


9.4(-03) 


...4<T^5cr^l7r 3 l7r 3 

ll y ll y 


V 




6.95 






1.6(+0) 


..Aal5a 2 A7ril7rl2Trl 

LL Lj LL L> LL 


V 


i 3 s+ 


1.68 








..Aal5allTT£lTr* 

u y u y 


V 




2.82 


2.84 






..A<J 2 h(T 2 llTllTTl 

u g u g 


V 


3 3 S+ 


4.30 








..Aalhall^lTT^l 

LL y LL y LL 


V 


4 3 S+ 


4.65 








..AallTrilirftnl 

LL LL LI LL 


V 


5 3 S+ 


5.63 








...5a2l4l7r 3 27ri 

y ll y ll 


V 


i 3 n 9 


0.95 


1.09 


0.82 




..Aa 2 u ba)lit A u liTl 

u, y ll y 


V 


2 3 n 9 


4.26 








..Aa l u ha 2 A^ u \4 


V 


3 3 n 9 


5.31 








..A<jl$o 2 g litilKpKl 


V 


4 3 n 9 

y 


5.96 








..Aalba 2 l^l-Kpnl 

LL Lj LL Lj LL 


V 


i 3 n u 


1.19 


1.37 


0.93 


4.6(-3) 


... 4^5(7;? It^Itt 3 

u g u g 


V 


2 3 n„ 


4.00 






2.2(-4) 


..AalBallTrllrf 

ll y ll y 


V 


3 3 n M 


5.10 






1.6(-2) 


..Aal5all7r£l7r 2 g 27ri 


V 


4 3 n n 


5.70 






2.8(-5) 


..AalballirilTTpTri 


V 


5 3 n u 


5.86 






1.2(-1) 


..Aa 2 u ha l g l^p^ u 


V 


1 3 A M 


2.77 


2.82 






...4a25^l7r 3 lvr 3 


V 


2 3 A U 


4.25 


3.56 






..Ao 2 J>o 2 g \^ u \it\2i<l 


V 


1 3 A 9 


6.67 








..A<j 2 u $a 2 g litllit 2 g 2irl 


V 



a Ground state energy including Davidson correction: -152.228684 Hartree 
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Table 9 

Vertical excitation energies, oscillator strengths and dominant configurations for 
I— C4H at the ground state equilibrium geometry 





Energy (eV) 


/ 


el 


Dominant 


Type 


State 


This 


Ref. 


This 


Ref. 


Configuration 






work 


(70) 


work 


(70) 






l 2 £+ a 


0.00 


0.00 






...9o l ln 4 2ii 4 


V 


2 2 £+ 


5.11 


4.76 


1.0(-6) 


<7(-7) 


...9cr 1 l7T 4 27T 3 37r 1 


V 


3 2 £+ 


7.36 


6.62 


9.5(-4) 


1.8(-3) 


...9o- 1 l7T 3 27r 4 37r 1 


V 


4 2 £+ 


7.74 




3.5(-l) 




Jr 7V 


R 


5 2 S+ 


8.76 




5.5(-l) 




...9tT 1 l7r 4 27r 3 3pi 

JT 7T 


R 


i 2 n 


0.37 


0.44 


7.8(-4) 


8.4(-4) 


...9a 2 lvr 4 2vr 3 


V 


2 2 n 


3.59 


3.31 


7.8(-4) 


9.8(-4) 


...9o- 2 l7r 3 2vr 4 


V 


3 2 n 


5.31 


4.71 


3.4(-4) 


7.2(-4) 


...9o- 2 lvr 4 27r 3 37r 1 


V 


4 2 n 


7.15 


5.92 


2.2(-3) 


9.2(-4) 


...9o- 1 10o- 1 l7r 4 2vr 3 


R 


5 2 n 


7.85 


6.82 


7.0(-3) 


8.8(-5) 


...9a 1 10a 1 l7T 4 27r 3 


R 


6 2 n 


7 QD 


7 Rd 


1.4(-2) 




...9a 1 10a 1 l7T 4 27r 3 


R 


7 2 n 


8.41 




2.2(-3) 




...9cr 1 llcr 1 l7r 4 27r 3 


R 


i 2 s- 


6.03 








...9o- 1 lvr 4 2vr 3 37r 1 


V 


2 2 S~ 


6.48 








...9o- 1 l7r 4 27r 3 3vr 1 


V 


3 2 S- 


8.58 








...9o- 1 lvr 4 2vr 3 3ri 


R 


4 2 S" 


8.86 








...9o- 1 l7r 4 2vr 3 3ri 


R 


1 2 A 


5.76 


5.20 






...9o- 1 l7r 4 2vr 3 3vr 1 


V 


2 2 A 


5.97 


5.21 






...9o- 1 lvr 4 2vr 3 37r 1 


V 


3 2 A 


7.85 








...9o- 1 lvr 3 27r 3 3ri 


R 


4 2 A 


8.10 








...9o- 1 lvr 3 2vr 3 3ri 


R 



a Ground state energy including Davidson correction: -152.228684 Hartree 
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Table 10 

Vertical excitation energies, oscillator strengths and dominant configurations for 
I— C5H at the ground state equilibrium geometry 





Energy (eV) 




/ 


el 


Dominant 


Type 


State 


This 


Ref. 


This 


Ref. 


Configuration 






work 


(73) b 


work 


(73) 6 






i 2 n a 


0.00 


0.00 








...llff 2 l7r 4 27r 4 37T 1 


V 


2 2 n 


3.05 


3.21 


5.1( 


-4) 


l(-3) 


...lla 2 l7r 4 2vr 3 37r 2 


V 


3 2 n 


3.91 


3.99 


1.8( 


-3) 


4(-3) 


...llo- 2 l7r 4 2vr 3 37r 2 


V 


4 2 n 


4.18 


4.19 


5.7( 


-2) 


3(-2) 


...lla 2 l7r 4 27r 3 37r 2 


V 


5 2 n 


5.20 


5.10 


3.5( 


-3) 


5(-3) 


...llo- 2 l7r 3 27r 4 3vr 2 


V 


6 2 n 


6.09 


6.35 


5.2( 


-2) 


l(-3) 


...lla 2 lA 4 3^ 2 


V 


7 2 n 


6.13 


6.11 


3.0( 


-3) 


1.4(-1) 


...lla 2 l7r 3 27r 4 3vr 2 


V 


1 2 S+ 


3.21 


3.19 


4.6( 


-3) 


9(-3) 


...llo- 1 l7r 4 2vr 4 3vr 2 


V 


2 2 S+ 


5.07 


6.28 


3.4( 


-3) 


l(-3) 


...lla 2 l^ 4 2^ 4 3s 1 


R 


3 2 S+ 


5.73 


7.01 


1.0( 


-6) 


2(-3) 


...ll(j 2 l7r 4 27r 4 3pi 


R 


4 2 S+ 


6.69 




1.9( 


-3) 




...lla 2 l^ 4 2^ 4 4s 1 


R 


5 2 £+ 


7.09 




2.0( 


-6) 




...llo- 2 l7r 3 27r 3 3vr 2 


V 


i 2 s- 


2.83 


2.80 


6.0( 


-3) 


7(-3) 


...lla 1 l7r 4 27r 4 3vr 2 


V 


2 2 5T 


5.12 


5.25 


1.6( 


-3) 


l(-3) 


...llo- 1 l7r 4 27r 3 3vr 3 


V 


3 2 S- 


5.99 


6.22 


2.6( 


-3) 


l(-3) 


...llo- 1 l7r 4 27r 3 3vr 3 


V 


1 2 A 


2.71 


2.62 


4.8( 


-3) 


3(-3) 


...lla 1 l7r 4 27r 4 37r 2 


V 


2 2 A 


6.06 


5.92 


2.6( 


-3) 


l(-3) 


...llcj 1 l7r 4 27r 3 37r 3 


V 


3 2 A 


6.63 


6.05 


1.5( 


-3) 


2(-3) 


...llcr 2 12o- 1 l7r 4 2vr 4 


V 



a Ground state energy including Davidson correction: -189.983327 Hartree 
b Their cc-p-VTZ+SP triple £ results if available, double £ results otherwise. 
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Table 11 

Dissociation and ionization energies" with respect to ground-state equilibrium en- 
ergy (in eV) for various species 



Species 


Products 


A-E'diss 


Ref. 




Ref. 


i-c 3 


C 2 + C( 3 P) 


4.63 


(75) 


12.1 


OS) 


l-C 4 


C 3 + C( 3 P) 


4.71 


(M) 


10.7 


(ZD 


l-c 2 H 


C 2 + H 


4.90 


TW 6 


11.4 


TW 


I-C3H 


C 3 + H 


3.27 


TW 


8.6 


TW 


C-C3H 


C 3 + H 


4.29 


TW 


9.6 


TW 


I-C4H 


C 4 + H 


4.65 


TW 


9.6 


TW 


Z-C5H 


C 5 + H 


3.56 


TW 


7.4 


TW 


HC3H 


C 3 H + H 


3.1 


(M) 


8.96 


(58) 


c-C 3 H 2 


C 3 H + H 


4.37 


(M) 


9.15 


(28) 


I— C3H2 


C 3 H + H 


3.87 


(M) 


10.43 


(78) 



a Experimental data refer to the adiabatic ionization potentials; computed values 
in this work to vertical ionization potentials 

b Computed in this work, using the same basis set and procedure for the products; 
dissociation energies were obtained at 20 Bohr 
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Table 12 

Photodissociation rates (in s _1 ) for various molecules in the unshielded interstellar 
radiation field 



Species 


Rate 


l-C 3 


5.0( 


-9) 


l-C 4 


8.5( 


-9) 


/-C 2 H 


1.6( 


-9) 


Z-C 3 H 


1.8( 


-9) 


c-C 3 H 


l.K 


-9) 


I-C4H 


3.7( 


-9) 


l-C 5 B 


1.3( 


-9) 


HC 3 H 


2.2( 


-9) 


c-C 3 H 2 


1.4( 


-9) 


I— C3H2 


5.1( 


-9) 



29 



